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Abstract 

A two-dimensional (2D) assembly of noninteracting, temperature-dependent, composite- 
boson Cooper pairs (CPs) in chemical and thermal equilibrium with unpaired fermions is 
examined in a binary boson-fcrmion statistical model as the superconducting singularity 
temperature is approached from above. The model is derived from first principles for the 
BCS model interfermion interaction from three extrema of the system Helmholtz free energy 
(subject to constant pairable-fermion number) with respect to: a) the pairable-fermion distri- 
bution function; b) the number of excited (bosonic) CPs, i.e., with nonzero total momenta — 
usually ignored in BCS theory — and with the appropriate (linear, as opposed to quadratic) 
dispersion relation that arises from the Fermi sea; and c) the number of CPs with zero total 
momenta. Compared with the BCS theory condensate, higher singularity temperatures for 
the Bose-Einstein condensate are obtained in the binary boson-fermion mixture model which 
are in rough agreement with empirical critical temperatures for quasi-2D superconductors. 
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1 Introduction 



Recent experiments indicate that composite bosons in ultra-cold clouds of most alkali atoms 
do indeed Bose-Einstein (BE) condense. Since Cooper pairs (CPs) of fermions (electrons or 
holes) in a many-fermion system form composite bosons in the sense of coupling to integer 
angular momentum, it is natural to consider the possible BE condensation of such pairs. The 
belief that some such condensate is central to superconductivity is more than 50 years old 
|§, ^ ^, 0]. High- Tc, as well as some organic, superconductors |§] are quasi-two-dimensional 
(2D). Quasi-ID superconductors have also been found |§. BE condensation (BEC) is impossible 
in two or less space dimensions [10| for usual or "ordinary" bosons (i.e., with a quadratic energy- 
momentum, or dispersion, relation). It is however still possible to have BEC in all dimensions d > 
1 for non-interacting bosons if they obey a linear dispersion relation [11| — such as CPs moving 
in the Fermi sea. This possibility arises because the Hohenberg theorem |10|, which prohibits 
BEC in 2D, relies on an /-sum rule based on the quadratic dispersion relation appropriate to 
bosons [^] moving in a vacuum. Such a linear dispersion relation for the CPs in a binary 
boson-fermion mixture model was recently found ||l^ to be consistent, without any adjustable 
parameters, with the anomalous linear (quadratic) temperature-dependence above Tc in the 
resistivity of optimally-doped (overdoped) cuprates whether hole- or electron-doped. For the 
observed quadratic T-dependence in overdoped samples linear-dispersion CP charge carriers are 
essential. 

Although extensive studies in the BCS-Bose "crossover" problem in superconductivity have 



spanned |14] a period of over thirty years, we note that BEC is distinct from the standard (i.e., 
zero center-of-mass momentum CPs) BCS theory condensation where only that one bosonic 
state exists. 

In this paper it is shown that in addition BEC is still possible in 2D even if the number 
of composite bosons (pairs of fermions) in a binary boson-fermion mixture is not fixed — as 
chemical/thermal equilibrium renders it coupling- and temperature-dependent — as long as the 
total number of fermions is fixed. This gives rise to an interesting statistical-mechanics problem 
irrespective of the particular mechanism for pair formation, and may have a vital application for 
superconductivity as well as for (neutral-atom) superfluidity such as in liquid ^He p|, dilute 
mixtures of ^He in ^He [16|, or in trapped Fermi gases [17|. The statistical model dealt with here 
may be seamlessly linked to BCS theory, via the fermionic energy gap, when boson/unpaired- 



fermion interactions are included as, e.g., in Refs. [^] and ||T^. However, in these two papers 
the quadratic CP dispersion relation has been assumed. The quadratic form has recently been 
shown |]20[| to apply only in the zero-density or vacuum limit when the Fermi sea disappears. 

In Section 2 we recall a 2D gas of fermions at T = interacting via a constant pairing 
interaction in an annulus about the Fermi surface — viz., the BCS model interaction. The binding 
energy of a single pair near the Fermi surface (CP problem) decreases practically linearly with 
the center-of-mass momentum (CMM) of the pair for all values of the momentum below breakup. 
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the breakup momentum typically being only about four orders of magnitude smaller than the 
Fermi momentum. In Section 3 we discuss why the interacting many-fermion system can be 
treated as a set of independent CPs (i.e., composite bosons with fermion number two) mixed in 
with pairable fermions which are not bound into pairs, i.e., unpaired fermions. In Section 4 the 
more realistic scenario is considered of the BEC of these pairs, incorporating pair breakup beyond 
a certain CMM. Although the number of pairs is not fixed but rather strongly coupling- and 
temperature-dependent, BEC is still possible in 2D. A simple binary boson- fermion statistical 
model is introduced by constructing the Helmholtz free energy for an ideal mixture of pairable 
but unpaired fermions plus paired fermions (both zero and nonzero CMM pairs), all in chemical 
and thermal equilibrium. The latter results through extrema of the free energy in: a) the 
pairable fermion occupation probabilities; b) the excited boson numbers (nonzero CMM CPs); 
and c) the ground boson number (zero CMM pairs). In Section 5 the coupling- and temperature- 
dependence of the boson number is derived. In Section 6 the critical BEC singularity temperature 
is obtained first by ignoring the unpaired fermions in a pure boson-gas model and then exactly 
for the boson-fermion binary mixture model from a T-dependent dispersion relation derived 
and calculated numerically, and results compared with empirical data. Finally, Section 7 gives 
conclusions. 



2 Cooper-pair dispersion relation 

Consider a 2D system of N fermions of mass m confined in a square "pen" of area and 
interacting pairwise via the BCS model interaction 

[ otherwise, 

where k = ^(ki — k2) is the relative wavevector of the two particles (see Fig. 3 below); Vk^k' 
the 2D double Fourier integral of the underlying non-local interaction ^(r, r') in the relative 
coordinate r = ri — r2; /u(T) the ideal Fermi gas chemical potential which at T = becomes the 
Fermi energy Ep = h^kp/lm with kp the Fermi wavenumber; Hujd = h'^k'jj/2m the width of the 
annulus about the Fermi circle in which the pairing interaction is nonzero, with ujd being the 
Debye frequency. This model interaction mimics the net effect of an attractive electron-phonon 
interaction overwhelming the repulsive interfermion Coulomb repulsions, whenever V > 0. 

If TlK = h(ki + k2) is the center-of-mass momentum (CMM) of a pair, let Ek be its total 
energy (besides the CP rest-mass energy). The eigenvalue (CP |^) equation for a pair of 
fermions at T = immersed in a background of — 2 inert, spectator fermions within a (sharp) 
Fermi circular perimeter of radius kp is then 

^ e{ki - kp) e{k2 - kp) , . 

V 2ek-{EK-n'K^/Amy ^ ^ 
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where again = h^k'^/2m, 0{x) is the Heaviside unit step function, and the prime on the 
summation sign denotes the conditions 

fci = |k + < (4 + kiy/^ and k2 = \k- < (4 + k^^/^ (3) 

ensuring that the pair of fermions above the Fermi "surface" cease interacting beyond the annulus 
of energy thickness hcoD in accordance with (|l|), thereby restricting the summation over k for 
a given fixed K. Without these restrictions (§) would just be the Schrodinger equation in 
momentum space for the pair. Setting = 2Ep — A^, the pair is hound if A/^ > 0, and (0) 
becomes an eigenvalue equation for the (positive) pair binding energy IS.K- Our /S.k and Aq 
should not be confused with the BCS energy gap A(T). 

Let A = g{Ep)V > be a dimensionless coupling constant with g{Ep) the electronic density- 
of-states (for each spin) at the Fermi surface in the normal (i.e., interactionless) state, which in 
2D is constant 



g{e) = L^m/2TTh^ = g. 



(4) 



The Cooper equation (Q) for the unknown quantity Ak is analyzed in Ref. |22]. For zero CMM, 
K 



0, it becomes a single elementary integral, with the familiar |21| solution 

2hL0D 



An 



p2/A _ I 



(5) 



valid for all coupling A. For small K, it is not too difficult to extract |22] the asymptotic result 



-[1 + 7r^(l + VTT^)]fiVFK + 0{K^) ^ Ao - -HvfK + 0(i^2 



F- 



(6) 



For 



where u = Qjj/Tp, and vp is the Fermi velocity defined through Ep = h kp/2m 
weak coupling, A — > 0, this linear dispersion relation gives the 2D analog of the 3D result stated 



as far back as 1964 in Ref. |23|, p. 33 (see also, Ref. |2J], p. 336) but with the 2D coefficient 
2/7r of the last expression of (pt) replaced by 1/2. 



3 Justification of boson formalism 

These CP boson-like structures could be called "quasi-bosons" since their creation and anni- 
hilation operators are known not to obey the usual boson commutation relations [^], p. 38. 
However, they do obey the Bose-Einstein distribution since the energy Ex of the CP is given 
only by the total CMM, but is independent of the relative momentum k. Thus, the possible 
energy states for the pair are Ek as defined in (§) . The number of pairs A^k that can occupy such 
a state can take on indefinite values since there exist also indefinitely many relative momenta, 
namely 

A^K = E-^k,K = 0,l,2,.... (7) 

k 
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Here, A/k,K = 0, 1 is the number of pairs characterized by both k and K, and is the same number 
as that characterized by definite ki and k2, namely A/i^k = ^ki'^k2 = 0; 1 where nk- = 0, 1 is 
the occupation number for a single fermion, these remarks ah referring to singlet pairing. Much 
of all this has been known p^ ] at least since 1958, albeit in somewhat different language. 

This view of an actual Cooper pair should not be confused with, say, an Anderson [^] 
phonon-like collective excitation (or modes) with weak-coupling dispersion relation — in 2D [p^ ] 
given by {l/V2)hvFK in the long- wavelength limit, and which evolves into the plasmon when 
Coulomb repulsions between fermions are switched on. CPs here, like deuterons, carry fermion 
number two and as such are definite in number (although in the CP case this number is coupling- 
and temperature-dependent) and can thus undergo BEC. This is distinct from collective exci- 
tations which are indefinite in number. Park |28|, e.g., distinguishes between "permanent" and 
"ephemeral" bosons, the latter sometimes being referred to as "quasiparticles" to distinguish 
from the former "particles" . 

For Nb ordinary bosons of mass ms and energy = Ti^K"^ /2mB in any positive di- 
mension, d > 0, a temperature singularity Tc |29] appears in the number equation Nb = 
'}2'kW~^'^~^^^^^^'^ — at vanishing bosonic chemical potential jiB < when the number 
of K = bosons just ceases to be negligible upon cooling. It is given by 

,2 r -\2/d 



ns 



9<i/2(l) 



rriBkB 

with ub the boson particle density Nb/L'^, and gd/2{^) the usual Bose integrals 

1 



(8) 



T{a) 



dx 



oo / 



z ^e^ 



1 



E^-tC(cx), 



1=1 



(9) 



where T{a) is the gamma function and ("(o") the Riemann zeta function of order a. The last 
identification in (^) holds when a > 1 for which C((t) < oo, while the series 5cr(l) diverges for 
fj < 1, thus giving Tc = for d < 2. For d = 3 one has C(3/2) — 2.612 so that (^) becomes 
the familiar formula Tc — 2>.?>lh?ri'^^ /niBkB of "ordinary" BEC. On the other hand, for bosons 
with (positive) excitation energy ek = — Ax given approximately by the linear term in (^) 
for all K, the singularity that lead to (S) now yields pH], for weak coupling, 



a{d)hvF 
kB 



d+l 

vr 2 ns 

r(^)<7d(i) 



i/d 



(10) 



where [11| a{d) = 1, 2/t: and 1/2 for d = 1, 2 and 3, respectively. Note that now Tc > for 
all d > 1, which is precisely the dimensionality range of all known superconductors including 
the quasi-lD organo-metallic (Bechgaard) salts [^. This is not inconsistent with the Hohenberg 
theorem |10| that there is no broken symmetry, i.e., long-range order, in a Bose fluid for d = 
1 or 2, since this is based on an /-sum rule for bosons with a quadratic dispersion relation. 
Indeed, both (0) and (ITol) are special cases of of the more general expression [|31| for any space 
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dimensionality d > and any boson dispersion relation ek 
constant, given by 



C.K"^ with s > and C» 



sT{d/2) {2TTynB 



s/d 



_27r'^/2r(d/s)5,/.(l) ■ ^^^^ 
In what follows the number of bosons will be temperature-dependent and it is in conserving 
the fermion number that the singularity arises. As is the case for the pure boson gas, a linear 
rather than a quadratic dispersion relation will be needed to obtain BEC in 2D. This emerges in a 
statistical model of an ideal binary mixture of bosons (the CPs) and unpaired (both pairable and 
unpairable) fermions in chemical equilibrium for which thermal pair-breaking into unpaired 
pairable fermions is explicitly allowed. 



4 First-principles statistical model 

Under interaction (||) at any T the total number of fermions in 2D is = L'^kp/2TT = A'^i -|- 
and is just the number of non-interacting (i.e., unpairable) fermions A^i plus the number 
of pairable ones A''2. The unpairable fermions obey the usual Fermi-Dirac distribution with 
fermionic chemical potential /j,. On the other hand, the A''2 pairable fermions are simply those 
in the interaction shell of energy width hujE) so that 

since the density of electronic states (§) is constant and the remaining integral exact. At any 
interfermionic coupling and temperature these fermions form an ideal mixture of pairable but 
unpaired fermions plus CPs that are created near the single- fermion energy /i(T), with binding 
energy A/^(r) > and total energy 

Ek{T) = 2^l{T)-/1k{T). (13) 

This generalizes the T = equation Ek = 2Ep — Ak introduced below (^). 

The Helmholtz free energy F = E — TS, where E is the internal energy and S the entropy, 
for this binary "composite hoson/pairahle-hut-unpaired- fermion system" at temperatures T <Tc 



is then [|2| 



F2 = 2 deg{e)\n2{e)e + kBT[n2{e)lnn2{e)+{l-n2{e)}ln{l-n2{e)}]\ 



/ fi—huj]j 
+ [2fi{T) - Ao{T)]Nb,o{T) 



+ Y,\[2f,{T)-AK{T)]NB,KiT) 

K>0 ^ 

+kBT[NB,K{T)l^NB,K{T) - {1 + AfB,i^(T)}ln{l + A^B,;^(T)}] |. (14) 
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The integral term is the contribution from the unpaired fermions and runs over all levels in 
the energy shell where the BCS model interaction is nonzero, n2(e) being the average number 
of unpaired but pairable fermions with energy e; the prefactor two comes from the spin. The 
second term gives the free energy of the bosons with CMM K = since their entropy is negligible 
in the thermodynamic limit; here Ns^iT) is the number of (bosonic) CPs with zero CMM at 
temperature T. The summation term represents the free energy of the bosons with nonzero 
CMM, while Nb,k{T) is that with arbitrary nonzero CMM K, and the cutoff Kq is defined ^] 
by Ai^Q = 0. The free energy F2 is to be minimized subject to the constraint that the total 
number of pairable fermions is conserved. 

If N2o{T) is the number of pairable but unpaired fermions, the relevant number equation for 
the pairable (i.e., active) fermions is then 

N2 = iV2o(T) + 2[Nb,o{T) + Nb,o<k<Ko{T)] = iV2o(T) + 2iVB(r), (15) 

where NBfi<K<Ko{T) denotes the total number of "excited" bosonic pairs (namely with CMM 
such that ^ < K < Kq), i.e., NBfl<K<Ko{T) = Y.o<k<Ko Nb,k{T). Minimizing the free energy, 
subject to the constraint that (15) be a constant, is equivalent to minimizing the grand potential 



^2 = F2-^X2N2. (16) 

a) Minimizing with respect to the fermion occupation probabilities n2{e) yields 
the Fermi-Dirac distribution with fermion chemical potential not fi, namely 

= + 1 ; p^{kBTr\ (17) 

Thus the total number of pairable (but unpaired) fermions then becomes 



N2oiT)^2 deg{e)n2{e) = 2 de , (18) 



deg{e)n2{e) = 2 ^ 

and should be compared with ( |l2|) for which contains only fi. Since in 2D g{e) is a constant 
(^, (|l8|) becomes the exact expression 



A^2o(r) = ^ln 



I _l_ ^-|3(^l~^l2~'huJD) 



(19) 



1 -|- g-l3{fj.-i^2+huir>) 

b) Minimizing 1^2 with respect to the excited boson numbers Nb,k{T), K > yields 
the Bose-Einstein distribution summed over all < ii' < Kq, namely 

Nb,o<k<Ko{T) ^ E Nb,k{T) = E [e/^^^^^^)-^'^^^ - 1]"^ (20) 

K>0 K>0 

The factor multiplying f3 in (pO|) may be rewritten as £k{T) — ^b{T), where £k{T) = Ao(r) — 
^k{T) > is a (nonnegative) excitation energy as suggested by (^, while plb{T) turns out to 
be 
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liBiT) = 2 [^2(T) - ^l{T)] + Ao(r). (21) 
This allows rewriting (|2^) in the more meaningful boson form 

Nb,o<k<Ko{T) = Y: [e^^^-(^)~^-(^)> - ir' (22) 

K>0 

where ^ b (T) is clearly the bosonic chemical potential associated with the entire binary mixture. 

c) Finally, minimizing Q2 with respect to the number of zero CMM (or, "ground 
state") bosons NBfi{T) gives 

2[^2(T) - fi{T)] + Ao(T) =0 (0 < T < T,), (23) 

valid only in the stated temperature range as NBfi{T) is negligible for all T > Tc. However, 
in view of (^Tj) this implies that Hb{T) = for all < T < Tc — which is precisely the BEC 
condition for a pure boson gas, even though one now deals with a binary boson-fermion mixture. 



5 Boson number 



To determine Nb{T) from (15) we need (^) which with (^) reduces to 



A^2o(T) = ^In 



1 + g-/3{Ao(T)/2-ft<^i3} 
1 + g-/3{Ao(r)/2+;ia;o} 



(0 < T < Tc 



(24) 



At T = two distinct coupling regimes emerge by inspecting (|2J): a) for Aq/2 < hujD or, from 
d) for A < 2/ln2 ~ 2.89, we have that A^2o(0) = 2g{^j){nujD - Ao/2); while b) for Ao/2 >nujD 
(or A > 2.89) A'^2o(0) is identically zero. Hence, the number of bosons Nb{0) at T = from ( p]5| ) 
is just Nb{0) = ^[^2 — -/V2o(0)]. Using (12) for A'^2 the fractional number of pairable fermions 
that are actually paired at T = 0, namely 2Nb{0)/N2 = 1 — N2o{0)/N2, becomes simply 



Ao/2hiOD = (e^/^ - 1)-^ — > e^2/A /f^^, A < 2/ln2 ~ 2.89) 
2Nb{0)/N2 = { ' V - / ; 

1 (for A > 2/ In 2 ~ 2.89). 

This fraction is plotted against coupling A in Fig. 1 as 2ns(0)/n2, since nB{T) = Nb{T)/L?' 
and 712 = N2/L'^. Since Nb{^) = ^fi'Ao for A < 2.89, only those fermions in an energy shell of 
width Ao/2 around the Fermi surface actually pair at T = 0, while for A > 2.89 all pairable 
fermions actually pair up since then Nb{0) = ghun = ^N2. This result contrasts sharply with 



the "heuristic model" |31], Eq. (|22[), where 2Nb{0)/N2 = 1 for all coupling, and is more in line 



with BCS theory which implies, in any d, a coupling-dependent fraction estimated (Ref. H] p. 
128; see also [||]) to be [g{EF)2A/2g{EF)nu;D? = {A/nuJof = (sinhl/A)-^ — > Ae''^'^, where 

A — *0 



8 



A = Uujd/ sinh(l/A) (again, not to be confused with the CP binding energy Aq) is the T = 
BCS energy gap for the same BCS model interaction (|l]) used in this paper; this is graphed as 



the thin curve in Fig. 1 and is seen to be much larger than (25) for fixed A. The breakdown of 
BCS theory for BCS model interaction couplings larger than A ~ 1.13 is clear both because: a) 
the alluded fraction cannot exceed unity; and b) physically, if the fermionic energy gap A > hujo 
no pairable fermions are available at all. This breakdown is indicated by the dashed curve in 
Fig. 1. (A strong-coupling many-body model differing from that of BCS theory but based on 
the BCS model interaction has been solved by Thouless p^). 

Also displayed in Fig. 1 are two finite-temperature results for 2Nb{T)/N2 = 1 — N2q{T) / N2 
which are obtainable from ( p^ ) for any T provided one knows Ao(T) for any T > 0. For T > 0, 
the 9{ki — kp) = ^(efci — Ep) in (Q) becomes 1 — n(^fcj, where n(^fc^) = (e^^*i + l)"-*^ is the 
Fermi-Dirac distribution with = ^fci ~ /^(^)) with the ideal fermion gas chemical potential 
//(T) in 2D being given exactly by 

Mr) = r'ln(e^^--l)^^i?i.. (26) 

Note that /u(r) decreases monotonically with temperature from its maximum value of Ep but 
does not turn negative until T = Tp/ In 2 ~ lAATp so that the BCS model interaction (|^), 
which requires ^{T) to be nonnegative, will not break down (i.e., become meaningless) over the 
entire range of temperatures relevant in this paper, see Fig. 4 below. Similar arguments hold 
for 6{k2 — kp). Since ki = k2 implies that = ^^3, (^) then leads to a simple generalization to 
finite-temperature of the K = CP equation, namely 

1 = A / d^e-"^ + 1)~2[2^ + Ao{T)]-\ (27) 





Its numerical solution for Aq{T) is illustrated in Fig. 2 for A = gV = 1/2 and v = huJp)/Ep = 
0.05. Note that if one assumes a T* such that Aq{T*) = 0, the resulting integral in (p7|) diverges 
and the equation can only be satisfied for A = 0; thus, there is no temperature T* at which 
"depairing" will occur for any fixed i' and any nonzero A. 

6 Critical temperature 

Neglecting the background unpaired fermions and modeling our system as a pure boson gas 
of CPs but with temperature-dependent number density nsiT), one converts the explicit Tc- 
formula (0) into an implicit one by allowing to be T-dependent. For d = 2 (^) becomes, 
since 52(1) = C(2) = tt^G, 

aVs hvp 



^^ = ^^\/--(^^)- (28) 
This requires ub (T) = Nb{T)/L? which in turn requires d24), along with Ao(T) as determined 
from d^), and is given by the expression 2Nb{T)/N2 = 1 - N2o{T)/N2. Solving (|2|) self- 
consistently with A = 1/2 gives the remarkably constant value Tc/Tp ~ 0.004 over the entire 
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range of = hcoD/Ep values 0.03 — 0.07 typical |Q of cuprate superconductors. On the other 
hand, the BCS formula T^^'^'^ ~ LlSezje-i/^ with A = 1/2 gives TjTp = 0.005, 0.008 and 
0.011 for v = 0.03, 0.05 and 0.07, respectively. Clearly, both sets of predictions are somewhat 
small compared with empirical cuprate values of T^/Tp that range from 0.01 — 0.1. 

To obtain the exact critical temperature without neglecting the background unpaired fermions, 
one needs the exact CP excitation energy dispersion relation exiT) = Ao(T) — Ak{T) which 
is neither exactly linear in K nor independent of T. To determine Ak{T) we need a working 
equation that generalizes Ref. ||2^ for T > via the new CP eigenvalue equation (|2^). Because 
of symmetry, see Fig. 3, one can restrict the angle 9 to the interval (0, 7r/2) where ki > k2, i.e., 



to quadrant I. Recalling (13), in d-dimensions (|2|) becomes 



\2ttJ J h'^{k^ -kl)/m + AKiT) + h^K^/4:m' ^ ' 

Here k^ is such that /u = fi'^k'^/2m and becomes fcp as T — > 0, while kf) is such that Hujd = 
Ti?k'^/2m. The prime on the integral sign now denotes the restrictions 

A:| = |k-iK|2 = k^ -kK cos 9 + \k'^ > kl, (30) 
A;f = |k+iK|2 = k"^ + kKcos9 + \K^ <kl + kl). (31) 

In Fig. 3 the darkest shading corresponds to these (BCS model interaction) restrictions. The 



conditions ( pO[) and ( pl|) can be studied separately but must be satisfied simultaneously. If 



K < 2^kl - A:!,, (|0D and (|3l|) are equivalent to 

{kl - sin2 0)1/2 + cos9 < k < [{kf^ + k^) - {K"^ sin^ 0] V2 _ cos 9 . (32) 



Note that for K > \/ k'^ + k'j^ — ^ k'^ — kj^ there exists a minimum value ^min of 9 given by 

K^/2{2kl + kl)-K^ 



COS 9^ir, = , ^ = , (33) 



while 0min = for < wA;2 _|_ /^2^ _ /^2 _ ^ introduce the dimensionless variables 



2(4 + A;2,)i/2 ^-fc^' Ep' "^-Tp-kl 



-)) ' r^F J^F '^^ 

with ksOo = fioJD = fi?k'j~)/2m and ksTp = Ep, where kB is Boltzmann's constant. Recall 
the d = 2 constant expression (^) for ^(e), the restrictions (^), and that for A' > and T > 
the step functions in @) 9{ki^2 - kp) = 9{\\¥. ± k| - kp) become [exp{-/3[7i^(iK ± k)V2m - 
/u(T)]} + \]~^ — but with 2efc in @) replaced by e^^ + e^j, Ep by /i(T) and Ak by Ak{T). One 
finally arrives at a working equation for the binding energy A^(T) that generalizes Eq. (18) of 
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Ref. namely 

1 = 1a dO Z"^™"^^^^^ + e^P{-/3K^ + (1 + 1^)'^^ + KicosO- 1]}]-1 



^ ^e_„ 2^2 + 2(1 + z.)k2 - 2 + A«(T) 

X [1 + exp{-/3[^2 ^ ^ j,)^2 _ 2^i + l> cos 6* - 1]}]-^ , (35) 



where U = htUD/fJ; £,min{0) = + ^ KCOs6 + a/1 - (1 + 2^)^2 8111^6', Cmax{0) 



VI + K COS 61+ J{1 + K^sin^e) and 



0min=l ° 2k< 1- V(l-Z^)/(l + l^), 

""^ \ cos"i(z^/{4Vl + z^ Kv/1 + 1^/2 - (1 + J^)k2}) otherwise. 
In (5^) we have introduced the more general dimensionless quantities ^ = k/k^, A^iT) = 



AK{T)/fi, where T = ksT/i^ oi $ = ^/5, and k = K/2^kl + k}). 

To obtain the critical temperature from the finite-temperature dispersion relation, besides 
solving (H) for Ak{T), one needs (0), (jl|), (|2|) and (|2|). At T = Tc both NsfiiTc) ~ and 
fJ-siTc) — so that one gets the implicit Tc-equation for the binary mixture gas 



l = ^ln 

V 



1 + g-{Ao(f,)/2-,.}/r, 



8(1 + z^) r'^oifc) K 



1 + g-{Ao(rc)/2+y}/Tc 

This must be solved numerically for the exact Tc for each A and v in conjunction with (^) for 
Ao(T') and ( |35| ) for both Ak(T') and k;o(7c)- Results for A = 1/2 are shown in Table 1 and Fig. 
4 for a range of values typical |Q of cuprates. 

In order to make this comparison we have taken T^/Tp — 1, a very good approximation up to 
the highest temperatures dealt with. For example, from Fig. 4 the highest T^jTp ~ 0.14 already 
gives T^/Tp ~ 0.9999 from (p6|), while for smaller T^jTp the values of T^/Tp are even closer to 
1. The Tc resulting from the exact dispersion relation for T = (dot-dashed curve) is somewhat 
higher than the exact result (full curve) but lower than that using the linear approximation for 
Ak{T) (dotted curve). It is also clear that the effect of using the exact or linear (in K) cases 
dominates the effect of the dispersion relation T-dependence. For cuprates d ~ 2.03 has been 



suggested [37| to be more realistic as it reflects inter-CuO-layer couplings but our results in that 
case would be very similar to those reported here for d = 2. 

Thus, for u = 0.05 the exact Tc is seen to be about 46% lower than the heuristic result 
found in Ref. |3|, Eqs. (|l|) and (H). It is curious that all results depend very weakly on the 
T-dependence of the CP binding energy Ak iT) , in spite of its being substantial throughout the 
temperatures spanned in this paper, as seen in Fig. 2. 

We defer study of the condensate fraction NBfi{T)/NB{T) below Tc and merely surmise that 
it may ultimately help explain the apparent absence [^8[ ^ in cuprates of the Hebel-Slichter 
peak of nuclear-spin (NMR) relaxation rates vs temperature for < T < Tc- Such a peak, 
originally seen (HO] in aluminum, is perhaps the most stringent and qualitatively convincing 
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experimental test of BCS theory (Refs. |23|, p. 71 and [41|, p. 79 ff). Besides cuprates, it 
is also absent |42| in several quasi-lD Bechgaard |^] and in several quasi-2D (ET) organic salt 
superconductors . 



7 Conclusions 

A simple statistical model treating CPs as non-interacting bosons in thermal and chemical equi- 
librium with unpaired fermions is proposed. The model gives rise to a boson number that 
is strongly coupling- and temperature-dependent. Since the CP dispersion relation is approx- 
imately linear, it exhibits a Bose-Einstein condensation of zero-CMM pairs at precisely two 
dimensions. Exact transition temperatures for the boson-fermion mixture based upon the exact 
CP dispersion relation are in reasonable agreement with empirical cuprate data. 

Needless to say, further corrections are yet to be included in the present simple binary 
mixture boson-fermion model, e.g., i) realistic Fermi surfaces, ii) Van Hove singularities |^] 
or other means of accounting for periodic-crystalline effects, as well as iii) the all important 
d-wave interfermionic interaction, iv) the boson-fermion interaction and v) residual interbosonic 



interactions. As to the latter, also generally neglected in BCS theory, if the lowering |44| of Tc 
in liquid "^He by about 29% with respect to the ideal Bose gas BEC is any guide, interbosonic 
interactions will also lower in a more realistic picture. As to the boson-fermion interaction, it 
is precisely this ingredient that enabled T.D. Lee and coworkers jl8], and Tolmachev |19] more 



generally, to link BCS and BEC through a relation stating that the BE condensate fraction is 
proportional to the (BCS-like) fermionic gap A(r) squared. 
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Table 1: Critical temperatures Tc/T^ for A = 1/2 depicted in Fig. 4 according to (|3^). The 
exact result is compared with the linear- in- approximation for both IS.k{T) and A;^(0) in 
order to test sensitivity of a temperature-dependence of the CP binding energy for nonzero K. 



V 


linear approx. with Ax(T) 


linear approx. with A/^(0) 


Exact 


0.03 


0.078 


0.068 


0.065 


0.04 


0.089 


0.079 


0.075 


0.05 


0.100 


0.088 


0.084 


0.06 


0.109 


0.096 


0.091 


0.07 


0.117 


0.104 


0.098 



Figure Captions 

Figure 1. Fractional number of pairable fermions that are actually paired, at three different 
temperatures, vs. coupling A for the present first-principles model (^) (thick curves) and 
estimated for BCS theory at T = as explained below ( p5| ) (thin curve). The number of 
pairable fermions with the BCS model interaction used is just ([T^); all of them are actually 
paired at T = in the heuristic EEC model, Ref. H Eq. (23). 

Figure 2. Temperature dependence oi K = ^ CP binding energy Ao(T) obtained numeri- 
cally from ( p7[ ) for A = 1/2 and v = 0.05. Note that when T = oo ( p7| ) is analytical for Ao(cxd); 
the latter then turns out to be about 10~^, so that the curve saturates from above to this value 
at T = cxD. 

Figure 3. Cross-section of overlap "volume" in momentum space (darkest shading) where 
the tip of the relative wavevector k (for two fermions with wavevectors ki and k2) must point 
for the attractive BCS model interaction (||) between them to be nonzero and form a Cooper 
pair of CMM magnitude hK. 

Figure 4. Critical BEC temperature Tc in units of Tp, resulting for the boson- fermion 
mixture from ( p6| ) for A = 1/2 for varying v = Tloj£i/^ ~ Qd/Tp: with no approximations (full 
curve); using Ak{T) evaluated at T = (dot-dashed); using the linear-in-i^ approximation for 
Ak{T) (dotted). The dashed straight line is the BCS formula Tc ~ 1.13606^/^^ for A = 1/2. 
The very lowest full horizontal line is the solution of the implicit Tc-equation (|2^) for the pure 
unbreakable-boson gas for i' = 0.03, 0.05 and 0.07. A huge Tc enhancement thus follows from 
the mere presence of background unpaired fermions. Cuprate data are taken from Ref. [5^]. 
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